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Abstract 

In the framework of SOS models, the dynamics of isolated and pairs of surface 
steps of monoatomic height is studied, for step-edge diffusion and for evaporation 
kinetics, using Monte Carlo techniques. In particular, various interesting crossover 
phenomena are identified. Simulational results are compared, especially, to those of 
continuum theories and random walk descriptions. 
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1. Introduction 

In recent years, the dynamics of steps of monoatomic height on crystal sur- 
faces has attracted much interest, both experimentally and theoretically. Ex- 
perimentally, equilibrium step fluctuations of isolated steps as well as steps 
on vicinal surfaces have been studied extensively, following pioneering mea- 
surements of Au(llO) [1] and Cu(lln) [2] surfaces. The step fluctuations are 
quantified by the equilibrium function G(t)= ((h(i,t + t eq ) — h(i,t eq )) 2 )i )teq , 
where h(i, t) denotes the position (or displacement) h of the step at site i and 
time t; one averages over reference times t eq in equilibrium and over the step 
sites. Typically, the experimental data could be fitted rather well by a power 
law, G(t) oc t a , with the exponent a being between about 1/4 and about 
1/2 [1-4] (attention may be drawn to a discussion on possible ambiguities in 
determining the intrinsic step fluctuations when using scanning tunneling mi- 
croscopy [5]). 



Theoretically, distinct atomic mechanisms driving the step dynamics have 
been identified, leading, at long times, to power laws with a being, indeed, 
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1/4 and 1/2 in the limiting cases of step-edge diffusion and evaporation- 
condensation kinetics, respectively. The scenarios have been established in 
Langevin (or continuum) [6-8] and other phenomenological [9] descriptions as 
well as in simulations on SOS models [6,10,11] (for the experimental relevance 
of the SOS model see, e.g., Ref. 12). In addition to G(t), the corresponding 
non-equilibrium function w(t) has been analysed, with w(t)= {(h(i,t) 2 )i be- 
ing the step fluctuations of an initially straight step, say, h(i,t = 0) = 0. w(t) 
is believed to show the same characteristic power-law at late times as G(t). 
Note that in the case of evaporation-condensation kinetics, an exact solution 
for w(t) for the closely related discrete one-dimensional Gaussian model, ap- 
plicable to isolated steps, with arbitrary number of sites is available [13]. 

In this article, we shall consider step-edge diffusion and evaporation kinetics 
of isolated and pairs of steps in the framework of SOS models. 

For isolated steps, effects of the step length, boundary conditions (periodic 
and pinned boundary conditions) and temperature are studied systematically, 
using Monte Carlo techniques, computing both equilibrium, G(t), and non- 
equilibrium, w(t), step fluctuations. Nontrivial and novel features at early and 
intermediate times as well as effects due to the finite length of the steps will 
be emphasized. 

For pairs of steps, the role of entropic repulsion on both types of dynamics 
will be elucidated. Studying various step quantities and correlation functions, 
intriguing crossover effects are identified, and our simulational results will be 
discussed in the context of continuum theories and random walk descriptions. 

The article is organized accordingly. In the next chapter, we shall present our 
simulational findings on isolated steps, followed by the chapter on pairs of 
steps. Each chapter is subdivided in discussing first evaporation kinetics and 
then step-edge diffusion. A short summary will conclude the paper. 

2. Isolated steps 



A surface step of monoatomic height may be described by the one-dimensional 
SOS model, defined by [14] 



h(i,t) is the step position at step site i and time t; the sum runs over neigh- 



H = eY,\h(i,t)-h(j,t)\ 




(id) 
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bouring step sites i and j, j = i ± 1. For a step with L active sites, i. 
e. a step of length L, the ends of the step may be either pinned, for in- 
stance, h(0,t) = h(L + l,t) — for pinning at equal step positions, or they 
may be connected by periodic boundary conditions, h(0, t) = h(L, t) and 
h(L + l,t) = h(l,t) (we set the lattice constant equal to one, with % and 
h(i,t) being integers). The time is measured in Monte Carlo attempts (MCA) 
or Monte Carlo steps per site (MCS), with 1 MCS = L MCA. 

A pair of steps may be described by two SOS models with the step positions 
hi(i,t) and fi2(i,t), where the subscript 1 refers to the, say, left step and the 
subscript 2 to the right one. To avoid a double step or crossing of the two 
steps, one assumes h 2 (i,t) > hi(i,t), leading to 'entropic repulsion' [15]. 

Step fluctuations, quantified by the equilibrium function G(t) or the non- 
equilibrium function w(t), result from the detachment or attachment of an 
atom at site i, decreasing or increasing the step position by one. We shall 
consider two types of dynamics: (i) a step atom is moved to a neighbouring 
site, i.e. a detachment at site % is followed by an attachment at site i± 1 (step- 
edge diffusion, s-d), and (ii) attachment and detachment events at the step 
are uncorrelated (evaporation-condensation kinetics, e-c). The probability of 
accepting the elementary move may be given by the Boltzmann factor of the 
energy change needed to execute that move, see Eq. (1), implying Glauber 
kinetics in the case of e-c and Kawasaki kinetics in the case of s-d [16]. Of 
course, other types of dynamics may be imagined, for instance, terrace diffu- 
sion [6,8,11], but they are outside the scope of this study. 

The step fluctuations may then be easily computed in Monte Carlo simula- 
tions. To calculate the non-equilibrium function w(t), one averages over an 
ensemble of N realizations (using different sets of random numbers), starting 
each time from the initial step configuration, say, h(i, t = 0) = 0. In deter- 
mining G(t), one first has to equilibrate the step, choosing a sufficiently large 
reference time, say, r . One may then generate a set of successive reference 
configurations, h(i, r m ), with r m = r + m5t, m — 0, 1, 2, ...M, 5t being a con- 
stant. G(t) = ((h(i,T m ) — h(i,t + T m )) 2 )i >Tm is obtained from averaging over M 
reference times r m and over the step sites i; it can be determined in a single, 
long Monte Carlo run ('dynamic averaging' [11]). When taking into account 
only one fixed reference equilibrium configuration, h(i, r ), and averaging over 
an ensemble of realizations (similar to the computation of w(t)), correlations 
will usually depend strongly on the choice of the reference configuration. 

To obtain accurate and reliable Monte Carlo data, rather extensive sampling 
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and care in choosing a suitable random number generator (for instance, the 
linear congruential random number generator may lead to erroneous results) 
are needed. 

In the following of this chapter, we shall summarize our main simulational find- 
ings on isolated steps, with comparisons to results on related exactly solved 
continuum or discrete cases [7,13]. Typically, in the simulations steps with up 
to 128 sites have been studied. 

At very early times, step fluctuations are diffusive [7,13], due to excitations 
at independent, equivalent step sites, i.e. w (or G) oc t. The diffusion coef- 
ficients in the non-equilibrium case are D™ ec = exp(— 2e//csT) for e-c, and 
D™ sd = 2exp(—4e/k B T) for s-d. The corresponding diffusion coefficients for 
the equilibrium function G are somewhat larger, with the enhancement factor 
tending to increase when lowering the temperature. Strictly speaking, the per- 
fectly diffusive behaviour holds, L > 1, only in the first Monte Carlo attempt. 
Already the second attempt leads to either sub diffusive or superdiffusive fluc- 
tuations. Accordingly, the effective exponent, e. g. for w, 

a&(t) = \n(w(t i+1 )/w(U))/ \n(t t+1 /U) (2) 

where t = y/t i+ iti, is then either smaller or larger than one. A simple calcula- 
tion shows that w becomes superdiffusive in the e-c case (in the second Monte 
Carlo attempt, there are local step moves, which cost no energy, at sites next 
to the initial step excitation), while it becomes subdiffusive in the s-d case 
(where the step moves in the second attempt are on average energetically more 
costly than in the first move). In contrast, the equilibrium step fluctuations 
G(t) tend to be subdiffusive in both cases. 

(i) For evaporation-condensation kinetics, w ec continues to be superdiffusive 
at early times before crossing over to the subdiffusive regime (see Fig. 1), 
where the step fluctations are governed by the line tension, approaching, at 
sufficiently long times for sufficiently long steps the form w ec cc t 1 / 2 (with the 
crossover time increasing like L 2 ) [6,7,13]. Similar properties hold for G ec . 

For finite steps, applying pinned (pd) boundary conditions, w ec will eventually 
reach the saturation value, = w ec (t — > oo), with the effective exponent 
cx™ s (t) decreasing rather sharply towards zero when approaching saturation at 
a time being again proportional to L 2 (see also Fig. 2 and Ref. 13). For suffi- 
ciently long steps, we find that w^ pd {L) approaches the continuum expression 
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Fig. 1. Effective exponent a™ s of the non-equilibrium fluctuations of an isolated 
step, showing the superdiffusive behaviour at early times, when using evaporation 
kinetics. A step of length L = 128, at various temperatures, with periodic boundary- 
conditions, was simulated, averaging over 10 7 realizations. 

[7] 

with the step stiffness £ = 2/c#T sinh 2 (e/2/cBT) [17,18] (here, attention may 
be drawn to the recent discussion on the experimental determination of the 
step stiffness [19]). For short steps, the length dependence of the saturation 
value is observed to deviate somewhat from the above linear L-dependence, 
due to lattice effects. Note that G e £ pd (L) = 2w™ pd (L). 

In marked contrast, when applying periodic boundary conditions, the time 
regime where the growth of the step fluctuations is governed by the line ten- 
sion, w ec oc t 1 / 2 , is followed by a regime in which the step fluctuations are 
dominated by the diffusive motion of the roughened step. The entire step now 
acts like a random walker, with the diffusion coefficient approaching a/L 
for large L, where a increases with temperature (similarly for G). The crossover 
time to the diffusive behaviour scales like L 2 [13].- Obviously, at T = oo, each 
Monte Carlo attempt will be successful, and one always encounters the initial 
diffusive behaviour of independently and randomly moving single step sites, 
Di. 
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(ii) For step-edge diffusion, applying pinned or periodic (pbc) boundary con- 
ditions, the step ffuctuations are bounded because of the conservation of the 
average step position in each move. The saturation value of w sd is closely 
related to that in the e-c case for pinned boundary conditions, see Eq. (3). 
We find w^ pd {h/2)w s ^ pd m (4/5)iu^ 6c for large L. For both boundary 
conditions, pd and pbc, G S ^(L) = 2w^(L). The time needed to approach the 
saturation value scales like L 4 . The asymptotic behaviour for indefinitely ex- 
tended, L — > oo, steps, w,G oc t 1//4 [6-10], is approached closely for finite 
steps at intermediate times. At earlier times, using periodic boundary condi- 
tions, the effective exponent displays an additional interesting non-monotonic 
behaviour [10j, both for w sd and G sd . As L increases, the corresponding max- 
imum in a™£ shifts to larger times (the shift being, possibly, proportional 
to L 4 ), and it gets weaker and weaker, approaching 1/4 from above. The 
maximum occurs in G sd about two times later than in w sd . For pinned bound- 
ary conditions, the corresponding effective exponents, for w sd and G sd , decay 
monotonically in time. 

3. Pairs of steps 

For pairs of steps, we studied both types of kinetics, evaporation-condensation 
and step-diffusion, applying especially periodic boundary conditions. 

(i) In the e-c case, we assumed the two steps of length L to be initially straight, 
hi(2)(i,t = 0) = h^L, and separated by d lattice units, i.e. — hf^ = d . 

In the simulations, the step length, 1 < L < 128, the initial separation dis- 
tance, 2 < d < 60, and the temperature were varied. 

Typical features of the dynamics of the steps are the early superdiflusive be- 
haviour, the step meandering and step roughening driven by line tension, the 
saturation of fluctuations due to the finite length of the steps, the diffusive 
motion of each of the roughened steps, and, in addition, phenomena result- 
ing from collisions between the wandering steps. To monitor these aspects, we 
recorded the non-equilibrium fluctuations of each step, w\ c (t) = wl c (t) =w{t), 
the (squared) width of each step b% c (t), k — 1,2, 

b e k c (t) = ((h k (i ,t)-h k (i c ,t)) 2 ) (4) 



where io refers to the first, io — 1, or last, io = L, step site, and i c refers to 
the center site of the step, with b\ c (t) = bl c (t) = b(t), as well as the average 
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distance d(t) between the two steps 

d(t) = j'£(h 2 (i,t)-h 1 (i,t)) (5) 

The brackets () in the definitions for b(t) and d(t) denote ensemble averages, 
i.e. averages over N Monte Carlo runs with different random numbers. 

By choosing suitable values for the step length L, the initial separation dis- 
tance do, and the temperature ksT/e, typical features of the step dynamics 
may be disentangled clearly, as exemplified in Fig. 2, showing the time de- 
pendence of the effective exponent a e ff of w, b, and d. There, with L = 128, 
d = 60, and ksT/e = 1.0, one can identify four characteristic successive times, 
see Fig. 2: ti denotes the time after which step fluctuations approach closely 
the simple law w(t) « c(t + t c i) 1 / 2 (that law is consistent with the asymptotics 
for indefinitely long isolated steps w oc t 1 ^ 2 ; most of the deviation of a™ s from 
1/2, as depicted in Fig. 2, may be attributed to the constant t cl ). At time t 2 , 
the width b(t) of each step begins to saturate, as reflected by the pronounced 
decrease in the corresponding effective exponent a b eS . The diffusive motion of 
the entire, roughened step, its width now being fully saturated, starts to dom- 
inate the step fluctuations, w(t), at time t 3 ; then w(t) ~ D^t + t c2 (in that 
purely diffusive regime at times t > t 3 the effective exponent of w may still be 
significantly smaller than one because of the constant t c2 ). The first collision 
between the two wandering steps is indicated by t 4 . The collisions, resulting 
in the entropic repulsion between the steps, give rise to the increase in the 
average separation distance d(t), and, at the same time, they lead to a slow- 
ing down of the step fluctuations, as measured by w. As for two non-crossing 
usual random walkers in one dimension (corresponding to the case L — 1), 
perfectly diffusive behaviour, with the diffusion coefficient D 2 of the isolated 
steps, shows up again asymptotically at later times. Thence, af s displays a 
dip after the first collision time (t 4 ), see Fig. 2, both for pairs of steps and two 
random walkers, L — 1. 

Of course, for other step parameters the characteristic crossover times may 
be less clearly separated, and/or more difficult to identify. For instance, when 
reducing the initial distance between the steps, do, or increasing the length of 
the steps, collisions between the steps may occur already much earlier than 
in the example discussed above. In that case, the separation distance d(t) is 
expected to grow like t 1 / 4 due to the entropic repulsion during the roughening 
of the two steps [20]. We simulated such situations as well. In any event, at 
late stages, one eventually encounters, for steps of finite length, the diffusive 
motion of the two steps acting like two non-crossing random walkers, as we 
confirmed in simulations. 



7 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 




- B— -n-e Bi 



10° 10 1 



l 2 | t 3 



10 10 
t (MCS) 



d * 



°o o ^ 



1 4 1 5 



10° 



Fig. 2. Effective exponent a e s for a pair of steps with periodic boundary condi- 
tions, using evaporation kinetics, of the non-equilibrium fluctuations w (circles), 
the width b (diamonds), and the distance between the steps d (squares). Steps of 
length L = 128, with do = 60, at ksT/e = 1.0, averaging over 18,000 realizations, 
were simulated. The lines are guides to the eye. 



(ii) For step-edge diffusion, the average position of each step is conserved. 
Equilibrium correlation functions are much suitable to describe the dynamics 
of the pair of steps [21]. In particular, we computed, doing dynamic averaging, 

C M (x,t) = ({h k {i + x,r m + t) - 4 0) )(^(^ r m) - h< i 0) ))i,r m (6) 



with k and I denoting either the same step (the 'intra-step-correlation func- 
tion' Cn = C 2 2 = C(x,t)) or different steps (the 'inter-step-correlation func- 
tion' C\2 = C 2 i = C s (x, t)), where i and i + x are step sites. Averages were 
taken over several reference times r m in equilibrium. 



Intra- and intercorrelation functions have been studied before in the frame- 
work of a continuum (or Langevin) description, dealing mainly with terrace 
diffusion [21]. In the s-d case, C s has been argued to vanish; this suggestion 
may be viewed with care, as indicated by our simulational findings for the 
discrete SOS model. 
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Fig. 3. Equilibrium intra-step-correlations C(x,t = 0), for isolated (full symbols) 
and a pair of steps with separation distance do = 5 (open squares) and 15 (open 
diamonds), using step-edge diffusion. The full circles refer to the continuum expres- 
sion, eq. (7). Steps of length L = 36, at ksT/e = 3.0 were simulated, using dynamic 
averaging over, at least, 10 6 reference times r m . 

In our simulations, we varied the step length, 10 < L < 48, the separation 
distance, 2 < d < 30, and temperature 0.8 < ksT/e < 3.0. In addition, the 
intracorrelation function C for isolated steps was computed. 

At t — 0, intracorrelations, C(x,t = 0), for an isolated step were observed to 
be close to the continuum expression [21] 

C(x,t = 0) = ^c(x/L) (7) 
with c(x/L) = 1/6 - x/L + (x/L) 2 , see Fig. 3. 

The rather slight deviations between the simulational data and the contin- 
uum description may be mainly attributed to the L-dependence in front of 
the scaling function c; the deviations diminish for longer steps. 

In the presence of the second step, we observe an overall reduction of the in- 
tracorrelations compared to those of an isolated step as long as the separation 
distance d Q is not much larger than the width, measured by b s = \Jb{t = oo), 
of each step. The reduction increases monotonically with decreasing d . In the 
limit of a pair of close-by steps, d = 2, the correlation function approximates 
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the form C(x, 0) oc Lc(x/L), reflecting the synchronization of the fluctuations 
of the two steps. However, for larger separation of the steps, the shape of 
C(x, 0) does not scale perfectly with x/L, albeit it does not deviate much, see 
Fig. 3, from the parabolic shape of the isolated step, Eq. (7). Indeed, we do 
not expect a scaling behaviour, because, at fixed temperature and separation 
distance do, b s and thence the 'overlap' of the fluctuations of the two steps 
depend clearly on the step length L. 

The form of the intercorrelation function at t — 0, C s (x, t — 0), is found to be, 
at small separation distance do, close to that of corresponding intracorrelation 
function C(x, 0), reflecting again the synchronization of the fluctuations of the 
two steps. The shape of C s changes towards a more (co-) sinusoidal form as d 
increases. Note that the correlations C s tend to increase with d , as long as 
the separation distance is comparable to the width of steps, b s . They decrease 
drastically when further enlargening do, because then the two steps only rarely 
encounter each other. 

As time proceeds, both the intra- and intercorrelations are described by, as 
may be obtained, e.g., from Fig. 4, 

C (s ) (x, t) = A cos(2nx/L) exp(-yt) (8) 

where the prefactor A and the exponent 7 depend on the type of correlations 
(intra or inter) and the specifics of the steps (temperature, step length, and 
separation distance). For instance, we find, for C and C s , 7 oc l/L 4 , with 
a rather weak (if at all) dependence on do- The exponential decay, with the 
observed L-dependence of the exponent, and the sinusoidal form may be ex- 
plained by applying Mullins' theory on the flattening of surface corrugations 
[22] to the steps. Indeed, the validity of that theory in describing the healing 
of perturbations in one-dimensional SOS models with Kawasaki dynamics has 
been checked before using Monte Carlo techniques [23]. 

According to Mullins' theory, for steps of length L with periodic boundary 
conditions, each step configuration, as it occurs at reference time r m , will 
eventually decay exponentially in time, with the relaxation time increasing 
like L 4 (in the s-d case), and it will take on a sinusoidal shape with the wave- 
length L [22]. The correlation functions, C(x,t) and C s {x,t), indeed, show just 
this behaviour. The time needed to reach the asymptotics depends, especially, 
on the content of other harmonics in the reference configuration. 

4. Summary 
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Fig. 4. Inter-step-correlations C s (x,t) for a pair of steps, using step-edge diffusion, 
with L = 24 and at k B T/t = 3.0, averaging over 200,000 reference times T m . The 
different curves refer, from top to bottom, to x = 0, 1,2, ...,L/2. Note that Eq. (8) 
is approached closely after about 10,000 Monte Carlo steps (MCS). 

Using Monte Carlo techniques, the dynamics of isolated steps and pairs of steps 
of monoatomic height has been studied for two types of kinetics, evaporation- 
condensation as well as step-edge diffusion, in the framework of SOS models. 
Periodic and pinned boundary conditions are applied. 



For isolated steps, in addition to the asymptotic behaviour at late stages, as 
obtained from Langevin (or continuum) descriptions, we observe interesting 
phenomena at short times, including the superdiffusive step fluctuations in 
the case of evaporation kinetics and a non-monotonic time dependence in the 
effective exponent in the case of step-edge diffusion with periodic boundary 
conditions. 



For pairs of steps, in the evaporation case, various time scales are identi- 
fied where distinct processes govern the step dynamics. At late times, ap- 
plying periodic boundary conditions, the two roughened steps act like two 
non-crossing random walkers. In the case of step-edge diffusion, results of a 
recent continuum theory have been checked and extended. In particular, the 
time-dependent correlation functions along each step and between the steps 
are non-zero, with their long-time behaviour being described by the classical 
theory of Mullins on the flattening of surface corrugations above the roughen- 
ing transition temperature. 
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